arXiv:hep-ph/0503251v3 26 Apr 2005 


The energy dependence of the saturation scale 

in DIS at low x* 


Masaaki Kuroda** 

Laboratoire de Physique Mathematique et Theorique, CNRS UMR 5825 
Universite Montpellier II, Montpellier 34095, France 

and 

Dieter Schildknecht 

Fakultat fiir Physik, Universitat Bielefeld 
D-33501 Bielefeld, Germany 
and 

Ludwig-Maximilians-Universitat Munchen, Sektion Physik 
D-80333 Munchen, Germany 


Abstract 

Consistency of the previously suggested color-dipole representation of deep- 
inelastic scattering (DIS) and vector-meson production at low x with DGLAP 
evolution allows one to predict the exponent of the W 2 dependence of the sat¬ 
uration scale, A 2 at (IU 2 ) ~ (W 2 )° 2 . One finds = 0.27 in agreement with 

the model-independent analysis of the experimental data from HERA on deep- 
inelastic electron scattering. 
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The present paper is concerned with deep inelastic electron scattering (DIS) 
at low x = Q 2 /W 2 <C 1. In short, we analyse the consistency between DGLAP 
evolution of the nucleon structure function F 2 (x, Q 2 ) and the color-dipole picture. 
We find that the exponent C 2 lcovy that in our formulation of the color-dipole 
approach determines the energy dependence of the total photoabsorption cross 
section at large Q 2 , a 7 * p (W 2 , Q 2 ) ~ ( W 2 ) C2 /Q 2 , or, equivalently, the energy 
dependence of the “saturation scale” A 2 at (PF 2 ) ~ (W 2 )° 2 , coincides with the 
result of previous fits to the experimental data, C 2 heory = (^® x P eriment _ 

For x = Q 2 jW 2 <C 1, the photon-proton interaction is dominated by the 
interaction of the photon with the quark-antiquark sea in the proton. The proton 
structure function for only contains the flavor-singlet quark distributions, 

and their evolution in Q 2 for Q 2 > Qq .TJ is in good approximation determined 
by the gluon structure function alone 2J, 


£>n(f,e 2 ) 

9 In Q 2 


Re+e~ 

97T 


a s (Q 2 )xg(x,Q 2 ). 


( 1 ) 


The notation in m is the standard one, we only note R e + e - = 3 J2fQ} — 10/3, 
where Qf denotes the quark charge and / runs over the contributing (rip = 4) fla¬ 
vors. In the physical picture underlying 0 the perturbatively calculable photon- 
gluon-scattering amplitude is supplemented by the gluon structure function that 
parameterizes the unknown proton properties to be determined experimentally. 
The structure function F 2 (x, Q 2 ) is related to the flavor-blind, flavor-singlet quark 
distribution 

xTj(x, Q 2 ) = n F (xq(x, Q 2 ) + xq(x, Q 2 )) ( 2 ) 

via 

F 2 (x, Q 2 ) = ^j^xE(x, Q 2 ). (3) 

In the color-dipole picture J], valid at low x 1 and any Q 2 > 0, in terms 
of the imaginary part of the virtual-photon-proton forward-scattering amplitude, 
the process of scattering proceeds via the fluctuation of the photon into a qq 
pair that subsequently scatters on the proton via (the generic structure of) two- 
gluon exchange J3J. The properties of the proton are contained in the color-dipole 
cross section 


r±,w 2 ) = J d 2 l ± a Mr (Q,W 2 )( (4) 


that depends on the two-dimensional transverse quark-antiquark separation, f±, 
[HI IS]- The function cr(l 2 , W 2 ) is associated with the gluon-transverse-momentum 
distribution in the proton, and the factor (1— exp(— il±r±)) in ([2]) is characteristic 
of the (QCD) gauge-theory structure. This factor originates from the couplings 
of the two gluons to either the same quark (antiquark) or to a quark and an anti¬ 
quark. Motivated by the mass dispersion relation of generalized vector dominance 
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u □ or, equivalently, life-time arguments JHj on the hadronic (quark-antiquark) 
fluctuation of the virtual photon, the energy W appears as the second variable be¬ 
sides r j_ or 1 1 _ in ( 0 . We refer to the literature Eiia for the explicit representation 
of the total virtual photoabsorption cross section in terms of the virtual-photon 
wave function describing the qq fluctuations of the photon and the dipole cross 
section in 0. 

At sufficiently large Q 1 2 , the dipole cross section in the limit of small interquark 
transverse separation, r 2 —► 0, becomes relevant. From 0 , for r 2 —► 0, 1 

<7(qq)p(rx, W 2 ) = Jf^TT J dQQa^t 2 , W 2 ). (5) 

By reformulating the 7 *-gluon-scattering approach underlying 0 in terms of the 
transverse position-space variable r±, one finds um that the gluon structure 
function xg(x,Q 2 ) in 0 is proportional to the right-hand side of 0, 


a,(Q 2 )xg(x,Q 2 ) = — j dl 2 T 2 a (qf)p (l 2 ,W 2 ) 


( 6 ) 


The gluon structure function, according to e, is proportional to the first moment 
of the gluon transverse-momentum distribution. 2 The validity of 0 is restricted 
to sufficiently large Q 2 , where both the notion of the gluon-structure function 
appearing in 0 as well as the f 2 —> 0 expansion of the dipole cross section in 
© are applicable. 

Actually only transverse and longitudinal ( qq)T = L (vector) states contribute 
to the imaginary part of the virtual forward-scattering Compton amplitude. The 
structure function F 2 (x, Q 2 ) may be represented EH in terms of the J = 1 pro¬ 
jections of the color-dipole cross section 0- The leading contribution to 


F- 2 (x,Q 2 ) = Ahff,, p (W'' 2 ,Q 2 ) = j^KjpfW 22 , Q 2 ) + <T 7l p(lV 2 , Q 2 )) (7) 


at sufficiently large Q 2 becomes 3 


F 2 (x,Q 2 ) = 


Q 2 

367r 2 


Fp+p— 


dl 


12 


4 1 


12 


4 Q 2 + ll 


/2 




1 Note that the scale for rq 2 that determines the validity of the expansion in 0 depends 
on the behavior of (t(l J 2 , W 2 ). Essentially, it is given by the effective or average value of If 
determined by W 2 ). 

2 Note that the Q 2 dependence of the gluon structure function at fixed x is contained in 
W 2 = Q 2 lx. This is at variance with the conventional assumption, where x occurs on the 
right-hand side in and the Q 2 dependence is intro6uced via the upper limit, Q 2 , of the 
integral in 0 . 

3 The expression (8) for F 2 (x,Q 2 ) is obtained from (3.13) and (3.14) in ref.[0 by expansion 
in powers of If 2 /(4Q 2 + If 2 ). 
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cr, 
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\j it a t n '> mi -4' 1 ,w 2 )), 


( 8 ) 


where l ' 2 is related to the gluon transverse momentum and the light-cone variable 
z via 

f , 2 


l ,2 = 


n c- (9) 

z(l - z ) 

Moreover, also the gluon structure function © may be represented in terms of 
the longitudinal part of the J = 1 projection of the color-dipole cross section, 

EU, 


a s (Q 2 )xg(x,Q 2 ) = 
In terms of the “saturation scale” 


8n 


dll 2 ll 2 a (q ^(ll 2 ,W 2 ). 


KAw 2 ) = 


7r 




a (oo) I “"-L "J- "(99) 


( 10 ) 

( 11 ) 


where the constant cd°°) will be explicitly defined below (compare (HUjl h (fTTHl 
becomes 

a,(Q 2 )xg(x,Q 2 ) = J-^a^AlJW 2 ). (12) 

So far our considerations have exclusively been based on the two-gluon ex¬ 
change structure embodied in the form of the color-dipole cross section ©. To 
proceed, we assume the flavor-singlet distribution B and the gluon distribution 
m to have identical dependence on the kinematic variables x and Q 2 . In our 
case, x and Q 2 appear in the combination W 2 = Q 2 /x. Both x£(x, Q 2 ) and 
a s (Q 2 )xg(x,Q 2 ) must then be proportional to A 2 at (IT 2 ). Since F 2 (x,Q 2 ) is pro¬ 
portional to x£(x,Q 2 ), compare ©, also F 2 (x,Q 2 ) in © must be proportional 

t°AL(W /2 )- 

Since, moreover, the longitudinal term on the right-hand side in © is propor¬ 
tional to the gluon structure in m and also the transverse contribution to 
F 2 (x,Q 2 ) in © must be proportional to A 2 at (kb 2 ). In terms of the integrals in 
©, the above requirement on the flavor singlet quark and the gluon distribution 
thus becomes 


12 7 /2 


dl±l_L(j( q q)J= 1 (l±iW) —r dl^l 


a (99-)f 


AA,w 2 ). 


( 13 ) 


Note that the constant r is related to the longitudinal to transverse ratio, 


^ lP (W 2 ,Q 2 ) 1 

^ P (W 2 ,Q 2 ) 2r 
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Our previous analysis [nm cairn of the experimental data on DIS was based 
on the equality of 




(15) 


i.e. on 


r — 1. (16) 

We found consistency with the experimental data, including m the available 
information on the longitudinal virtual photoabsorption cross section. We will 
henceforth put r = 1. Upon inserting (El) and intoducing A 2 at (PU 2 ) from (ITU) . 
F 2 (x,Q 2 ) from (JHJ) becomes 


F 2 (x,Q 2 ) = 


Rp 


367T 3 


■cr 


(0O) A IJW 1 


1 

1 + 2 


(17) 


where the sum on the right-hand side refers to the sum of the transverse and 
longitudinal parts. 

We now insert F 2 (x,Q 2 ) from (El) and the gluon distribution ED into the 
DGLAP-evolution equation E to find the interesting constraint 


or, alternatively, in terms of the observable F 2 (x, Q 2 ), 

—-—F 2 (2H/ 2 ) = -FMV*). 
dlnW 2 ' ’ 3 { ’ 


(18) 


(19) 


We adopt a power-law ansatz for A 2 at (lU 2 ), 


W 


c 2 


W 


2 \C 2 


Ai,(H' 2 ) = B [—* + 1 — B 7772 


w. 


W r 


( 20 ) 


identical in form to the one previously employed El El HU Ej in (successful) 
hts to the experimental data on the virtual photoabsorption cross section, 


a rp (W 2 ,Q 2 ) = a rp ( V (W 2 ,Q 2 )), 


( 21 ) 


and on deeply virtual Compton scattering (DVCS) in terms of the scaling variable 

dE| 4 

^.<3 2 ) = ^t. (22) 


Ait (W 2 )' 


4 For the sake of clarity, we introduced the notation Ag at (lh 2 ) = A 2 (W 2 ) for the quantity 
previously denoted by A 2 (W /2 ). 
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Inserting the power-law ansatz EQD into m , we deduce 

1 fl\ C2 

a - = 3(2) • (23 > 

and accordingly 5 

Cf eory = Q>276> (24) 

The theoretically deduced magnitude of the exponent = 0.276 in El 

is in agreement with the model-independent fit [nun to the experimental data, 

C 2 X p = 0.275 ± 0.06. (25) 

The model-independent fit assumes that <r 7 * p (W 2 , Q 2 ) may be represented by a 
smooth function of the scaling variable ?](W 2 ,Q 2 ) in (12211 with A 2 at (kb 2 ) from 
P|). As a consequence of the generality of the ansatz, any theoretical prejudice 
with respect to the empirical validity of scaling in r/(W 2 , Q 2 ) is excluded. 

The fit based on an explicit ansatz for the color-dipole cross section (compare 
El below) gave the more precise result 

Cf p = 0.27 ± 0.01 (26) 

in agreement with our theoretical result El- 

To summarize: the choice of W as the relevant variable in the color-dipole 
approach together with the requirement that the singlet quark distribution and 
the gluon distribution have identical dependence on the kinematic variables, 
W 2 = Q 2 /x in our case, converts the DGLAP-evolution equation ED into a 
constraint that allows one to predict the exponent C 2 . The agreement with 
experiment supports the validity, at least as a relevant approximation, of the 
underlying assumptions. 

The choice of W 2 = Q 2 /x as the relevant variable in our approach was mo¬ 
tivated by the color-dipole approach and its generalized-vector-dominance inter¬ 
pretation. We note that the dependence d 

F 2 ~ (W 2 f 2 (27) 

is closely related to the so-called singular solution of the gluon evolution equation, 
where 6 

F 2 ~ In Q 2 x~ x ~ ^(W 2 r = (^ 2 ) A (28) 

{Q ) 

with A > 0.25 being fixed and equal to the input value at all Q 2 . I 11 a restricted 
but relevant range of Q 2 , EZD is similar to ESD- 

5 The arguments leading to iTHli were implicitly used in the third paper of ref. m without, 
however, fully realizing their significance. 

6 Compare the discussion in section 7 of ref. m and the literature quoted there, in particular 
ref. (T7i. 
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The conventional application of the color-dipole approach {3J to DIS does not 
explicitly introduce the J = 1 projection of the color-dipole forward-scattering 
amplitude. One starts by an assumption on a(l 2 ,W 2 ) in ©, rather than its 
J = 1 projection, d( g g)j=i(/ ± /2 , IT 2 ). We briefly elaborate on how our approach 

can be formulated in terms of d(/j_ 2 , W 2 ). 

Consider the ansatz 

?(/?, W 2 ) = ^-6 (T 2 - |aL,(W’ 2 )) • (29) 

Inserting (ESI) into ED, we immediately recover the gluon structure function (Ill- 
Evaluation of the J = 1 parts of El yields 

■?(«)£-(C, W*) = h,T(l±, AL(M /2 ))«('? - l^JW 2 )), (30) 

the explicit form of the function _/ 2 , As at (IE 2 )) being irrelevant in the present 
context. We only note the normalization of (USD 




£ =1 (l± iW 2 ) — / dll a( q qy=i{ll ,W ) — 


a 


(oo) 


7r 


(31) 


Evaluating the longitudinal part of F 2 (x,Q 2 ) in © by substituting (© 3 ) , we 
recover our previous result m for the longitudinal contribution. The normal¬ 
ization EH) suggests to approximate El by 


cr 




qvMKiW 2 ) = 


a 


(oo) 




= — s(K 2 - a Uw 2 )). (32) 


7T 


With El inserted into (©, we now obtain not only the longitudinal, but also the 
transverse part of F 2 given in (ED- The direct evaluation of F 2 , inserting El, 
however, yields 


F 2 (x,Q 2 ) 


Re+e~ 

367T 3 




fli e ® 2 

12 n al,(w 2 ) 



(33) 


in distinction from ED- Our requirement of identical singlet quark and gluon 
distributions that is contained in ED, ED and El is not fulfilled by the ansatz 
El- It is not known whether an ansatz for cr(l 2 , W 2 ) can be given such that 
ED with r = 1 be valid. For the time being, we have to accept the equality 
as a valid approximation for the J = 1 projections that describes the experimen¬ 
tal data on DIS in the low x diffraction region with the predicted value of the 
exponent C 2 . 

Encouraged by the result E3D, we now examine the coupled system of equa¬ 
tions for singlet quark and gluon evolution. 
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Substituting the gluon structure function m and the singlet quark distribu¬ 
tion x11{x, Q 2 ) from (0) with (ED, into the evolution equations, with the power- 
law (El for A 2 at (kb 2 ) and the notation t = In Q 2 , one hnds 


dAL(W 2 ) 

dlnW 2 


A s 2 at (hh 2 ) 


dy 


2vr 


P qq(y) + P qg{y ) V 


c 2 


(34) 


and 


dA Lt( w2 ) 

din W 2 


= A L(w 2 ) 


1 da s (t ) 


+ KAW Z ) 


a s (t) 

2\ 


dt 

r 1 


2tt 


dy 


+ 


a 




7r 


Pggiy) + P gg{y) y 


c 2 


(35) 


The hrst equation, CT . without relying on the approximation contained in the 
right-hand side of (HD, describes the evolution of the flavor singlet quark dis¬ 
tribution, while the second equation, El, describes the evolution of the gluon 
distribution. By noting that <9A 2 at (lT 2 )/<91nl/b 2 = C 2 A 2 at (lT 2 ), and upon evalu¬ 
ating the integrals on the right hand side, in El, we obtain 


C 2 


0.044a s (t) + 


C 2 + 3C 2 + 4 

2(C 2 + 1)(C 2 + 2)(C 2 + 3) ’ 


(36) 


The numerical value of 0.044 in the (small) (^-dependent correction proportional 
to a. s {t) in El was obtained by inserting C 2 = 0.276. Solving El f° r C 2 , we 
find 


C 2 = 0.044a s (t) + 0.260 

^ 0.265, (37) 

upon disregarding the (weak) Q 2 dependence of a s at large Q 2 and inserting 
a s = 0.11. 

A similar approach, when applied to the gluon-evolution equation El , leads 
to 

C 2 = —— J + C' 3 (C 2 )a 2 (t) + C, i (C 2 )oi s (t) 

a s (t) dt 

^ 0.275 (38) 

The dependence of the coefficients C' 3 (C 2 ) and C 4 (C 2 ) on C 2 is directly calculated 
from El- The value of C 2 = 0.275 was obtained by consistently solving El , 
using a s (t) = 0.11. 

As a result of our analysis of the complete evolution equations, by comparing 
El and El with El , we conclude that the power-law ansatz El with a con¬ 
stant (Q 2 -independent) value of C 2 hcoiy of magnitude C 2 heoiy = 0.276 according 
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to (J2H) is consistent with evolution. The additional a s -dependent contributions 
in (EH which have been ignored in (JTJ) hardly affect the value of the exponent 

c 2 . 

So far in this paper we were concerned with DIS at low x 1 and sufficiently 
large Q 2 3> Ag at (hh 2 ) where the QCD-improved parton picture and the color- 
dipole picture are dual descriptions of the underlying physics. 

For details on the extension to Q 2 <C A 2 at (W 2 ) and the (successful) description 
of the experimental data based on the ansatz EH, we refer to refs. nsmunsi- 
We only note the Q 2 —> 0 limit in addition to the large-Q 2 limit, 


CV P (^ 2 ,Q 2 ) 


a r MW\Q 2 )) = 

In - 


a 


3n 


R P 


c-O 


(°o) 


Q 2 +m 2 - 


AL.(^ 2 ) 

2 Q 2 


(■Q 2 « AL(iv 2 )), 
('Q 2 » ALGv 2 )). 


(39) 


With A 2 at (W 2 ) from (ITU) , the asymptotic limit of Q 2 A 2 at (W 2 ), in (jHHll coin¬ 

cides with EH. We also note 


cd 00 ) = 7f 


W 2 ). 


(40) 


Hadronic unitarity requires a ( - 00 - ) to at most show a weak W dependence. The fit 
to the experimental data led to cd 00 ’ = constant. 

It is worth noting that the virtual photoabsorption cross section, or F 2 (x, Q 2 ), 
for x = Q 2 /W 2 -C 1 and any Q 2 only depends on the integrated quantities in 

EH an d (1401) . 

From EH, at any Q 2 , for sufficiently large energy, such that A 2 at (kF 2 ) 3> Q 2 
and r]{W 2 ,Q 2 ) 0, the virtual photoabsorption cross section approaches the 

“saturation limit” of PHUH 


lim 

VF 2 —>oo 
Q 2 fixed 


vrp(v(W 2 ,Q 2 )) 

a lp (W 2 ) 


1. 


(41) 


The quantity A 2 at (IF 2 ) indeed sets the scale for the limiting behavior (1411) . The 
terminology “saturation scale” for the effective gluon transverse momentum squared, 
(l/6)A 2 at (!F 2 ), originating from the underlying two-gluon exchange, is indeed ap¬ 
propriate. 

In conclusion: the previously formulated color-dipole approach to DIS (the 
generalized vector dominance/color dipole picture, GVD-CDP) has been exam¬ 
ined with respect to the underlying singlet quark and gluon distribution. Both of 
these distributions being proportional to the saturation scale, A 2 at (lF 2 ), we find 
that the evolution equations lead to a remarkable constraint on the value of the 
exponent of the W 2 dependence that agrees with the experimental result. 
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